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We present a systematic characterization of the radio frequency (RF) spectra of homogeneous, paired atomic 
Fermi gases at finite temperatures, T, in the presence of final state interactions. The spectra, consisting of 
possible bound states and positive as well as negative detuning {v) continua, satisfy exactly the zeroth- and 
first-moment sum rules at all T. We show how to detect the i^ < Q continuum arising from thermally excited 
quasiparticles, which has not yet been seen experimentally. We explain semi-quantitatively recent RF experi- 
ments on "bound-bound" transitions and, thereby, predict the associated effects of varying temperature. 
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The superfluid and normal phases in trapped Fermi gases 
undergoing BCS to BEC crossover are presenting us with 
novel forms of superfluidity. An important characteristic of 
the superfluid is the pairing gap which is best probed using 
radio frequency (RF) spectroscopy IjJ, |2|] . This technique has 
been applied experimentally in a trap integrated yll a and to- 
mographic |2] fashion. While early theoretical work fl |5|l 
addressed trap effects, more recently attention has been on fi- 
nal state effects L6, J,^^ J^ although, unfortunately, only at 
low or zero temperature. Many have viewed the importance of 
these experiments as a means of quantitatively measuring the 
ground state pairing gap, thereby testing different approaches 
to BCS-BEC crossover. Our point of view is that finite tem- 
perature is crucial to full experimental understanding as well 
as reliable assessments of theory. 

The goal of this paper is to present a single formalism for 
the RF spectra at all frequencies and all T, including final state 
effects. A successful theory of a Fermi gas near unitarity, not 
only (i) has a pairing gap which appears fllsl lltl at T* > Tc 
but which, as T is decreased, (ii) exhibits a second order phase 
transition, at Tc- Studies of this smoothly varying (from above 
T* to T = 0) pairing gap, reminiscent of its counterpart in the 
high Tc superconductors, may elucidate some of the physics 
of the cuprates [11]. On physical grounds |3, S 0] it is clear 
that the RF current I{i') reflects the pairing gap A(r) rather 
than coherent superfluid order. At odds with this observation 
is the fact that all crossover theories which include pairing 
fluctuations 11121 [l3l [l4l [l5ll except the present one, lead to 
first order transitions at Tc. In a related fashion, alternative 
calculations OHIgllial of I{i') consider only the low or zero T 
superfluid and/or separately the normal phase even though, at 
T < T*, the presence or absence of superfluid order in the RF 
spectra should not lead to fundamentally different physics. 

We consider a homogeneous system which is relevant to 
recent tomographic experiments |2] At T j^ 0, the spectrum 
consists of (possibly) bound state contributions which either 
appear at positive or negative detuning, v and, (always), pos- 
itive as well as negative ly continuum contributions which re- 
flect the pairing gap, and can be used to measure its size. We 
emphasize the i^ < continuum which derives from thermally 
excited quasiparticles has not yet been seen experimentally 
nor addressed theoretically. A central finding is that it can be 
strongly enhanced by final state interactions and made visi- 



ble in future tomographic experiments. Near unitarity, final 
state effects make it possible to extract (using sum rules) the 
gap A as well as the chemical potential /i. We explain semi- 
quantitatively recent low T experiments and make predictions 
for the accompanying temperature dependences which should 
be observable. 

The RF technique focuses on the three lowest energy 
atomic hyperfine states (two of which are involved in the pair- 
ing, while a third provides a final "excited" state for one com- 
ponent of a pair). For definiteness, we first consider a su- 
perfluid of pairs in the equally populated hyperfine 1-2 lev- 
els and apply a radio frequency LU23 to excite the atoms in 
state 2 to state 3, as described by a Hamiltonian given in 
Refs. nil S 0, II H- The RF response function can be ob- 
tained following the standard linear response theory |4, 5, 6]. 
Here we formulate the finite T, RF problem using a diagram- 
matic scheme which can be made compatible with the dia- 
grams in Ref. 8, although attention in that paper was restricted 
to very low temperatures. We will see below that our diagram- 
matic scheme reduces at T = to the approach of Ref. §1 
This correspondence, and indeed, all diagrammatic formula- 
tions 1 7] of the RF experiments are based on a T-matrix ap- 
proach. The T-matrix used here (for the 1-2 channel) is con- 
sistent lUullTJl with the BCS-Leggett ground state equations 
and involves one bare and one dressed Green's function. We 
have 

tl2(Q) = 9i2 + E Gi{K)Gl{Q - K) (1) 

K 

ti^Q) = 9ii + E Gi{K)Gl{Q - K) (2) 

K 

where we have introduced the dressed Green's function G = 

[{G°y^ - S]"^ and G" is the Green's function of the non- 
interacting system. Here the subscripts indicate the hyper- 
fine levels, K = {iuji, k), Q = (ifi„, q) are 4-momenta with 

J2k = ^X); X^k' ^'■''■' '^^'^ ^i ^^'^ ^« ^'"^ fermion and bo- 
son Matsubara frequencies, respectively. Throughout we take 
h ^ kB = I and assume a contact potential (so that the strict 
Hartree self-energy vanishes) and a (nearly) empty population 
in the hyperfine 3 state so that G3{K) « G'^{K). As has been 
demonstrated elsewhere lUlll . it is reasonable to take the self 
energy (on the real frequency axis) in the Green's functions 



Gi and G2 to be of the generalized BCS form 

A2 



E(w,k) 



Ck 



(3) 



although this approximation is not essential. Similarly, 
we have shown llllll that, below T^, A(T) is constrained 
by a BCS-like gap equation which can be written as 1 + 
512X12(0) = where xi2(g) = T.kGi{K)GI{Q - K), 
in conjunction with a fermion number equation. More gen- 
erally, the propagator for noncondensed pairs is of the form 

ti2(Q)=gi2/[l + .gi2Xi2(Q)]. 

We emphasize a distinction between the pairing gap (which 
we call A) and the order parameter, called Asc. The differ- 
ence between these two energy scales can be shown llllll to be 
associated with noncondensed pair effects parameterized by 
the pseudogap A^g defined by 



AlJT) = A\T)-^l,{T). 



(4) 



Here we note that A^^ — — Sq=^o *i2(Q)' which allows Tc 
to be determined llllll as the temperature where lS.sc first van- 
ishes. We find Tc — 0.25Tf at unitarity. It is convenient 
notationally to define a form of Gor'kov F function in terms 
of the pairing gap as 



AG2{K)G'i{-K) 



A 



El 



F{K), 



where Ek = VCfe + A2(T), and Cfe = ffc - A*. Cfe = k^/2m. 
Because of the constraints imposed by the BCS-like gap equa- 
tion, ti2{Q) diverges at Q = so that it is reasonable to set 
Q in ti2 to zero, i.e., ti2{Q) « — (A^/T)(5(Q). This assump- 
tion, which leads to the simple form of Eq. (O, is not essential 
for understanding the physics but it does greatly simpUfy the 
calculations flq]. 

The resulting diagram set for the RF response function, 
D{Q), is shown in Fig. 1 and this last approximation is equiv- 
alent to treating the Aslamazov-Larkin (AL) diagram (called 
Dal) in Fig- 1 at the BCS mean-field level, leading to the op- 
posite momenta ±K for particles 1 and 2 in the diagram. The 
leading order term, Do{Q), of the response function appears 
as the bubble on the left and was introduced in Ref. 4. The 
term on the right, Dal{Q), depends on A, not Asc, and in- 
corporates final-state effects via the interactions 312 between 
1 and 2 and 313 between 1 and 3. We neglect the effects aris- 
ing from the interaction between 2 and 3. This is consistent 
with the approach in Ref. 6. This second term has appeared 
previously in studies of the superfluid density lUTT . 

Writing out the AL diagram yields 

Dal{Q)^ [Y,FiK)Gl{K + Q)'\\i3{Q) . (5) 

K 

For the RF field, Q = (ifi„,0) so that D{inn) = D{Q). 
We take ns satisfying fi^k^a) = 0, where ^^,3 = ek - A*3- 
Then the RF current, given by the retarded response function. 



is /(zy) 
find 



-(l/7r)ImZ?^(J7), where il = i^ + zi — /i3, and we 



DiQ) = DoiQ) 



[D2{Q)r 



m/4TTai3 + Di{Q) ' 



(6) 
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Figure 1: (Color online) Feynman diagrams for the RF response 
function D(Q). The left bubble is the lowest order Do, whereas 
the right diagram, Dal, is associated with final state effects. Here 
thin (thick) lines stand for bare (full) fermion propagators, the dashed 
line for ii2, approximated as the condensate, and double wiggly line 
for ii3. The numbers in blue indicate the hyperfine levels. 



and t^^{Q) = m/47rai3 + Di{Q), where ai3 (and ai2) are 
the s-wave scattering length in the 1-3 (and 1-2) channels, 
respectively. Here DoiQ) = Ek G2{K)GI{K + Q) 



lfiEk)-fi^k,3) 



^ \-i^n+ Ek — ^k,3 i^n — Ek — S.k,3 



l-f{^ k,3)-f{Ek) 2 
i^k 



(7) 



and/o(i^) = -(l/7r)Imi:>^(r2). We also define DalQ) 
Y.kF{K)GI{K + Q) 



E 

K 



A 
2E'k 



l-I{Ek)-I{^k,3) f{Ek)~f{^k,3) 



. ifln — Ek — S,k,3 

Y.kGi{K)GI{Q 



-fit 



k.3) 



1 



iVLn+ Ek - Cfe,3 

K) - Ek(l/2e;.) = 

/(Cfe,3) - I[Ek) 



(8) 



v-^ in 



andi)i(Q) 

v^ \f{Ek) 

j^l inn-Ek-ik,3 "^' inn + Ek-^k^3 "^J ^ ^^ 

(9) 
After analytical continuation and change of variables, we have 
il ± E]^ — £,k.3 ~ f^ i £^k ~ ^/c- Importantly, the denominators 
here are the same as those which appear in ti2- Furthermore, 
at i^ = 0, f{^k.3) is cancelled out so that 

*r/(0) = {gi3 - 912) + tuiO) = 9ii - 9i2- (10) 

It follows that the complex functions Dq{Q), Di{Q), and 
Z?2(Q) are the same as their wave function calculation coun- 
terparts fv] when the pairing gap A is chosen to be order pa- 
rameter Asc and T = 0. It is i/ not 51 that should be identified 
with the experimental RF detuning. 

After some straightforward algebra, one can show that 
when (713 = gi2 there is an exact cancellation such 
that I{i') ^ ^ii^)- In general, we have Io{i') = 



1,B., 



(l/7r)(AVj^^)Imt~3''"(i/), and 



/M = 



1 
512 
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513 



4(^) 



itr^^'^^HP 



47ra 



13 



•-13 

m 

47rai2 



■A^ 



ImtfgM, (11) 



where if3(i/) =tfji{n). 

Equations (fTTT i are a central result of this paper which make 
it clear that final state effects in the RF current directly reflect 
the T-matrix in the 1-3 channel. In general, features in the RF 
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Figure 2: (Color online) RF current J(i/) as a function of RF detun- 
ing u for transitions from unitarity l/kFa.12 = at 834 G to final 
state (a) l/kpUi^ = —1 and (b) —0.5 in the BCS regime, corre- 
sponding to Tf —31 and 124 kHz, respectively. The temperatures 
are T/Tp =0.1 (Black solid), 0.3 (red dashed) and 0.45 (blue dot- 
dashed lines). Here Tc = 0.25Tf- The sharp lines next to the right 
continuum in (b) correspond to bound states. Inset: Lowest order RF 
current 1 0(1^) vs z/. 
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Figure 3: (Color online) RF current /(//) as a function of detuning u 
for 1 — > 2 transitions in a 1-3 superfluid of Tf — 40 kHz (a) from 
l/fcFflis = -0.804 to final states l/kFa2s = at 811 G, and (b) 
from l/kFai3 = -0.524 to l/fcFa23 = 0.68 at 750 G, for different 
temperatures as labeled. Here Tc/Tf ~ 0.15 and 0.17, respectively. 
In (b) when T is high and A is small, the two peaks around u = 
may not be resolvable experimentally. 



spectra derive from the poles and imaginary parts of Eqs. (|6]l- 
(|9]l. The spectrum may contain a bound state associated with 
poles at j/Q in ii3> as determined by t^-]'{i'o) = 0. This leads to 
the so called "bound-bound" transition. In addition, there is a 
continuum associated with both the numerator and denomina- 
tor in the first of Eqs. ( fTTT ). with each contribution spanned by 
the limits of i^ = ^fe ± Ek, i.e., -{^/JJFT~^ + ^i) <v <Q 
and V > y^fi^ + A^ — /i. The continuum at positive frequen- 
cies is primarily associated with breaking a pair and promot- 
ing the state 2 to state 3. This represents the so-called "bound- 
free" transition. On the negative detuning side, the continuum 
is primarily associated with promoting to state 3 an already 
existing thermally excited 2 particle. The spectral weight of 
the negative continuum vanishes exponentially at low T as 
g-A/T Therefore, there is a strong asymmetry in the con- 
tinuum with the bulk of the weight on the positive frequency 
side for low T. If the bound state falls within the negative 
continuum, it will acquire a finite life time, and decay quickly 
at high T. 

Of importance, in assessing a theoretical framework for 
computing the RF current are the two sum rules associated 
with the total integrated current and the first moment or "clock 
shift" 1 6]. Using the Kramers-Kronig relations between Re tf^ 
and Im ifg, it is easy to prove that, not only in the ground state, 
but also at finite temperature, Eq. ( fTTI ) satisfies 



both diagrammatic contributions are included. It is easy to 



dvl{v) 
dv V I{i') 



n2-n3, 



Att 



1 
ai2 



1 
ai3 



(12) 
(13) 



where 712 and n3(= 0) are the density of state 2 and 3 atoms, 
respectively. In this way we find for the clock shift 



A2 



J dvvl{v) 
J dv I{v) n2 — na 47r 



1 
ai2 



1 
ai3 



, (14) 



which agrees with Ref. |6] when 713 — > 0. It should be stressed 
that this sum rule is satisfied only when ai3 ^ and when 



show that at large v, Iq (i/) 



-3/2 



, Imtf; 



"1/2, so that 



I{i') ^ v~^l'^, in agreement with Ref. 8. Clearly, the first 
moment of liy) is integrable, whereas the first moment of 
/o (v) is not. Finally, Eq. ( fTTT i reveals that the spectral weight 
(including possible bound states) away from 1/ = will dis- 
appear when the gap A vanishes. 

Figures 2(a) and 2(b) illustrate the behavior of the spectrum 
I{y) when the initial state 1-2 pairing is at unitarity (i.e., at 
834 G) and the final state 1-3 pairing is on the BCS side of the 
1-3 resonance, for temperatures T jTp = 0.1, 0.3, and 0.45. 
The parameters we use are taken from Ref. [l^. The inset of 
Fig. 2(a) indicates the behavior in the absence of final state 
effects for the same temperatures. The asymmetry of the con- 
tinuum around v = Q, discussed earlier, is evident even in this 
leading order bubble diagram. As T is raised the spectrum be- 
comes more symmetric. In contrast to the findings in Ref..la, 
and as consistent with experiments Ij]] on trapped gases, we 
do not find a substantial pairing gap at T jTp « 1. In Fig. 2(a), 
the final state interaction Xlkpaiz = ^1 is relatively weak, 
and there is no bound state. In contrast, at l/fci?ai3 = —0.5 
(or Tp w 6/iK) in Fig. 2(b), a bound state emerges at low 
T (although it disappears at moderate temperatures when the 
gap becomes small). For the low Tp ^ 2.5/iK used in Ref.Hl 
we do not find a bound state. These results are consistent 
with r = calculations of Basu and Mueller |9]. It should be 
stressed that, at 834 G for a typical Tp, when the bound-bound 
transition occurs, it is barely separated from the asymmetric 
bound-free continuum, which is always present. 

Figure 3 presents the analogous plots at different T for RF 
transitions from an initial 1-3 superfluid with Tp — 40 kHz at 
(a) 811 and (b) 750 G, which are on the BCS side of the 1-3 
resonance (which appears at 690 G). The system is subject to 
an RF field promoting state 1 to state 2. "Bound-bound-like 
transitions" 11911 now appear In Fig. 3(a), the bound state falls 
within the negative detuning continuum. Importantly, the dis- 
appearance of the bound state with temperature is preceded by 
a very unusual two-peaked spectrum in the negative detuning 
regime, which is seen at the two higher T . We can under- 
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Figure 4: (Color online) RF current I{v) as a function of detuning 
u for a 1-3 superfluid with RF excitation from state 3 to state 
2. The black curves are calculated at experimental parameters of 
il/kFai3,l/kFai2,T/TF) = (0.4,3.3,0.2), (0.0,2.6,0.1), 
(-0.3,2.0,0.1), (-0.7,1.1,0.09), (-0.9,0.6,0.09), and 
(-1.2,0.0,0.06) from low to high fields. The red dashed 
curves are calculated at twice the temperatures. The sharp lines on 
the left indicate bound states. For comparison, experimental data are 
marked by arrows for bound peak locations and by triangles for the 
continuum. 12011. 



stand this unusual structure as a combination of the peak from 
the negative continuum which appears very close to ly ~ 0, 
(as also seen in Figure 2) and the near-by bound state peak. 
At even higher T, the spectral weight will shift almost com- 
pletely to the region near j/ = 0^ and the bound state decays 
rapidly. As j/ ^ 0^, the negative continuum peak is a com- 



bined effect of the vanishing Imt{| and the diverging factor 
1/j/^ in Eq. ( fTTT i. In Fig. 3(b), the bound state is outside the 
continuum, and the binding energy is fairly insensitive to tem- 
perature. We have chosen experimentally accessible parame- 
tesrs here, so that the unusual double-peaked structure in /(i^) 
at i^ < should be observable. Finally, we emphasize that the 
highest T cases in Figs. 2 and 3 are at or above Tc, so that the 
continuum appears only because there exists a pseudogap in 
the fermionic spectrum. 

Figure 4 addresses recent data [2] associated with 1-3 pair- 
ing and RF excitation from state 3 to state 2. The calcula- 
tions of I{iy) shown in the (black) solid curves in all six pan- 
els were performed with experimental parameters, and should 
be compared with Fig. 4 of Ref. ||31. To help in the compar- 
ison a number of data points (normalized to the same peak 
height) have been inserted. The sharp bound states will, in 
the data, be broadened both instrumentally and from limited 
spatial and energy resolution. Except for a slight broadening 
which we have ignored here, our calculated black solid curves, 
which incorporate final state effects, can be seen to be in semi- 
quantitative agreement with experiment. We anticipate that at 
higher T (red dashed lines), the negative ly continuum states 
should start to become apparent. Despite the presently good 
agreement, we feel the ultimate test of any theory must involve 
a test of its predictions, such as those shown here. 

At unitarity, the best way to measure A(r) is using the 
sum rule in Eq. ( fT4b and its experimental counterpart. To- 
gether with the t/ > continuum threshold which appears at 



VA^ 



M^ 



/i, one can also determine /i and hence the fac- 



tor (3. This analysis is possible only in the presence of final 
state effects. Because I{i') at general T depends on the total 
pairing gap A(r), the size of the order parameter and pseu- 
dogap cannot be separately inferred (except when the order 
parameter vanishes above Tc). 
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